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Abstract 

The studies of influence of spin on a photon motion in a Schwartzschild spacetime 
is continued. In the previous paper [2] the first order correction to the geodesic 
motion is reduced to a non-uniform linear ordinary differential equation and the 
equation obtained has been solved by the standard method of integration of the 
Green function. If each photon draws a world line specified by this solution then 
light rays from infinitely distant source form a caustic which does not appear without 
the spin-gravity interaction. The goal of the present work is to obtain explicit form 
of caustic. 

Keywords: Worldline of photon; Spin-gravitational interaction; Schwarzschild 
spacetime. 

1 Introduction 

In our recent works [TJ [2] deviation of photon world line from a null geodesic was ob- 
tained in the first order approximation on spin-gravity interaction in Schwarzschild 
space-time. If all world lines form a parallel beam from an infinitely distant source, 
spin- gravitational interaction of photon should give the following effect. Since the 
world lines obtained do not lie wholly on one semi-plane of constant azimuthal angle 
if they do not cross the axis of symmetry on which the source is, but pass at some 
non-zero distance from it. Then some neighbohood of the back semi-axis is shad- 
owed by the deviation. The shadow has certain boundary formed by the envelope 
of the rays, hence, this surface appears as caustic of this family of curves. The goal 
of the present work is to find out the shape of the surface. 

To obtain the shape of the surface we need to consider some details of calculations 
of the reference geodesic and our first-order approximation. Therefore it is useful 
to recall the techniques and the denotions used. We use the standard denotions: 



coordinates for the Schwartzschild space time are {t, r, 9, <p} and null geodesies lying 
wholly on the 6 = pi/2 equatorial "plane" are presented in the parametric form 

dR dR „ 
where 



'E 2 r 2 - L 2 (l - rJr) 

R(r) = s(t)l \ *' ' dr, (2) 

J ro r(l-r g /r) 

e(t) = signumi which follows from the corresponding solution ^ = Et — Lip + R(r) 
of the Hamilton-Jacobi equation < d*S>,d^ >= 0. The constants are chosen such 
that under t = each geodesic reaches the minimal value of r, r = tq. Then under 
t — > — oo the geodesic approaches the starting point at the source and under t — * +oo 

it runs away to the flat asymptotics. We put (p\ t _>_ 00 = 0, so that (fo = 

It must be noted that unlike action of a massive particle which can be used as a 
parameter on the geodesic, that of massless one ^ takes constant values on null 
geodesies in question. Therefore, ^ cannot be used as a parameter on the geodesies. 

The form of 0-geodesics determined by squaring (pQ) is given in terms of hyper- 
elliptic functions. In fact, we do not need the exact expressions because in the 
commonplace approximation we employ: r g <C D where D = L/E stands for the 
impact parameter, all expressions simplfy (see, for example, Ref[3 ). The radial part 
([2j) of the action function ^ can be represented as 



flW = - E(s) /„7(T^ 

Following [3] we substitute 



D 2 



o V (1 - r g /r) 2 r 2 (l - r g /r) 



dr. 



r(r — r g ) = r' 2 , so that r' ~ r — r g /2. 



Then the function R(s) becomes 

r-rg/2 



2r > m dr'. 



ro~rg/2 

Moreover, we have 

r' 2 + 2r g r' — D 2 = x 2 - B 2 ', B 2 = D 2 + r 2 , x = r' + r g « r + r g /2. 
Another substitution in the integrand yields: 

rr+r g ( % 2 _ 1)2)1/2 

R(s) » -e(s)E / ^ — dx. 

Jd x-r Q 



This result allows to represent the shape of light ray approximately in analytical 
form considered below 



2 Analytical representation of the shape of the 



ray 

We see that the factor annullating the integrand at the lower limit is separated and 
now we can safely expand another factor putting r g /x < r g /D <C 1 that gives: 

rr+r g l % 2 _ £)2\l/2 

R(s)^-e(s)E ^ '- — (l + rjx)dx. (3) 

JD X 

It is seen that the integral admits analytical representation. 

The shape of the null geodesic is determined by dependence of angle (p on s given 
by eq. ([I]) due to Hamilton-Jacobi theorem. Substituting eq. ([3]) into the eq. for (p 
we obtain: 

rr+r a /2 ( D r q D 1 

The constant parameter ifQ can be eliminated by corresponding choice of initial 
value of this coordinate: <p(— oo) = accepted in our previous work [2]. This choice 
yields the following explicit form of (p as function of the variable s: 



ip(s) « 7r/2 + r g /D + e{s) < arccos h 



D . r g ^/x 2 - D 2 



I x Dx 



x=r+r g /2 



Asymptotical shape of the caustic to be found is formed by s > parts of the 
geodesies: 



V(s > 0) « tt/2 + r g /D + arccos ° {r I 9,2) + — . ( i ) 



It was shown in our previous work [2] that in the first order approximation the 
geodesies deviate only in dg direction. The deviation is specified by the vector of 
the deviation 5x = 5x 2 r~ 1 de. Its only non-zero component for the outgoing branch 
of the trajectory (s > 0) has the following asymptotical behavior under r — > oo: 

^ - & (5) 

where we have introduced local Cartesian coordinate x 2 such that dx 2 = rdO. Fi- 
nally, approximate world line of photon is given by the reference geodesic and the 
small vector of deviation just provided. 



3 The shape of asymptotical part of the caus- 
tic 

Since the shape of caustic of stationary beam of light is stationary as well as the 
Schwarschild space-time itself, we do not need the time coordinate on the world line 
of given photon and can consider only the shape of separate light ray in the space 
endowed with standard spherical coordinates {r, 92, 8}. 

The envelope of a beam of such rays is just the caustic to be found. Besides, axial 
symmetry of the envelope allows to represent it by the shane of curve along which 



the surface crosses equatorial semi-plane < ip < vr/2, 9 = tt/2. In other words, to 
obtain the desired surface it suffices to find out the distance between envelope and 
the axis p 2 = r 2 sin 2 ^) + (5x 2 ) 2 as function of the coordinate r. This task reduces 
the problem of constructing the envelope of family of curves on the semi-plane 

p(r, D) = yjr 2 sm 2 ((p) + (6x 2 ) 2 , (6) 

where each curve is labeled with certain value of the impact parameter D. Value 
of ip is given by dU). It is convenient to use the coordinates {r, p} in the semi-plane 
instead of the standard ones {r, <p}. 

Appearance of caustic can be explained geometrically as follows. Due to the 
spin-gravitational interaction the rays instead of crossing the axis pass at distance 
p from it. The rays which are asymptotically straight lines, under given value of 
the impact parameter D lie on one-sheet hyperboloids whereas the geodesies do on 
cones with same axis. Unlike the cones, the hyperboloids of distinct D constitute 
a family which has an envelope. In two-dimensional picture they are hyperbolas 
whereas the geodesies are straight lines incident to the axis. The hyperbolas have 
an envelope we are constructing. In this section we simplify the task assuming that 
the envelope of hyperbolas almost coincides with the curve on which their points 
closest to the axis lie. We call this curve "simplified" caustic and first, we explore it. 
Afterwards we show that in our approximation this curve coincides with the genuine 
caustic constructed in classical approach. 

Since all the results are obtained under assumption that the value r g /D is small, 
the impact parameter is to be taken sufficiently big D 3> r g that allows to obtain only 
asymptotical behavior of the caustic. To do it we need explicit form of the expression 
for the function sin<p(r) in our approximation. Straightforward substitution of the 
equation (JH) yilds: 

sinfa) « D(l/r + r g /r 2 ) - 2r g /D « - - (7) 

It is convenient to represent the envelope by the curve r = r x (D), p = p(r x (D), D), 
where r x {D) is the value of the coordinate r under which the ray with impact 
parameter D crosses the axis. To find it we solve the equation 

sin <p(r, D) = 

which due to equation ([7]) reduces to 

D 2 

r* = a- (8) 

This allows to express the variable p in the neighborhood of the crossing point 
r = r x : 

where A is wavelength. In other words, asymptotically the envelope becomes a 
cylinder of radius which disappears in astrophysical scales. 



4 Classical approach 



Since simplified caustic is too narrow to be observed it is necessary to obtain the 
genuine caustic. The shape of the genuine caustic can be extracted from shape of 
envelope of curves ([6]) on which the following determinant vanishes: 



dp 
dD' 



Indeed, if we use the parameter D as a coordinate instead of p, each curve is given 
by D = const. Since the curves have an envelope the Jacobian of the transformation 
vanishes on it, consequently, the envelope to be found is zero of the determinant. 
Thus, the equation of the envelope can be written as follows: 



dp 


dp 




dp 


dp 


dr 
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dr 
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dr 


dD 





dD 



0. 



(10) 
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It is seen that under 8x 2 = no envelope appears and since 5x 2 
envelope can appear only at distances of order ([8]). Referring to the equations (|5I6|7[) 
we have: 

J^T 2 T 2 



r 2 r 2 



+ D Z + 



4r g r. 



p E 2 D 4 D 2 
Now, substituting the expression into the equation (|10p we obtain: 







dp 
dD 



+ 2D-8 



,^2 ~,2 
9 



D 3 D 3 ' 

Taking account that 1 / E ~ A <C D we have finally the same result as in the equation 



2^ 

r Rj 

4 r 2' 
9 

for r being the coordinate of point of intersection of the ray and the axis. 



5 Angular size of the caustic 

Usually, angular size of an object is specified by the angle under which the rays 
from it diverge on the axis. In this sense the caustic has zero angular size because 
its shape is cylindric. However, the rays are not parallel on it, consequently, some 
caustic must be seen. Below we calculate the angle under which they pass near the 
axis. Since the rays are straight near the crossing point the angle is 8(j> = 

ar r=r x 

Substituting this into the equation @ to the above expression gives: 



r 



dip _ D _ D _ 2r g _ [2r, 



' 9 



The angular size is doubled angle 5ip, namely, 2^J2r g /r where r is coordinate of 
observer. This result seems to be somewhat unexpected because angular size of 
caustic obtained this way does not depend on the wavelength whereas in the limit 
of zero wavelength the caustic does not appear at all. 



Since angular size of the caustic obtained this way is valid for all wavelengths it 
can be obtained also for geodesies which expose no spin gravitational interaction. 
The angle in question apparently if formed by geodesies focused onto the axis and 
must be observed because any telescope detects only direction from which radation 
comes. Naturally, radiation in this case comes from directions which consitute a 
cone which must be seen as a bright ring. However, this image has nothing to 
do with spin-gravitational effect and exposes only the fact of lensing. As for the 
spin-gravitational interaction, it does not produce any special observable effect but, 
probably, some diffusion of the bright ring with dependence of the diffusion on the 
wavelength. 

6 Conclusion 

This is the closing work of the series started with the article [4] in which Papapetrou 
equation was derived for a massless spining vector particle. In the work [T] derivation 
was slightly improved and an attempt was made to obtain an approximate solution 
of the equation. The final version of derivation of Papapetrou equation from La- 
grangian of electromagnetic field is presented in our work [5] . In our work [2] it was 
shown that the method used in the work |lj is valid only on the first half of the 
photon world line and another method, valid on the whole world line is proposed. 
The new method allows one to obtain the whole world line as a small deviation from 
the reference geodesic. 

All this work was completed with the purpose to explain why black holes appear 
in observations as bright rings. We hoped to prove that these rings are images of 
caustics formed by photons inclined from geodesies by the spin-gravitational inter- 
action. This would be another test of general relativity. In fact, it is not so. As was 
shown above, caustics produced by this interaction are asymptotically too narrow 
and have actually zero angular size, therefore, is not observable. Instead, there exists 
ring-shaped image of another nature which has nothing to do with spin-gravitational 
interaction. Angular size of this image does not depend on the wavelength and ex- 
poses only focusing of geodesies. As for the spin-gravitational interaction, it can 
only give a subtle spectroscopic effect on the image. 
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